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1 .  Introduction 


Radio  signals  that  pass  through  the  Earth’s  atmosphere  are  sensitive  to  electrically 
neutral  particles  that  cause  a  time  delay  in  the  signal  depending  on  the  dielectric  constant 
of  the  media  and  path  length.  Additionally,  if  the  signal  does  not  travel  in  parallel  with  the 
density  gradient  of  the  atmosphere,  it  will  be  subject  to  atmospheric  refraction,  causing  its 
path  to  deviate  from  a  straight  line.  Geometrical  bending  of  the  signal  path  is  greatly 
intensified  in  the  lower  atmosphere,  most  notably  the  troposphere,  due  to  the  curvature  of 
the  air  layer  near  the  Earth’s  surface.  This  bending  effect  on  atmosphere  is  significant  at 
low  (ground  station-to-spacecraft)  elevation  angles  (less  than  3  degrees).  Consequently, 
an  improved  modeling  technique  (exponential  tropospheric  model  or  Choi  model)  for  an 
excess  propagation  delay  due  to  the  atmosphere  must  be  incorporated  into  the  geolocation 
measurement  software  in  order  to  take  frill  advantage  of  propagation  delay  and  to  reduce 
the  effect  of  elevation  angle  dependent  errors  induced  by  projecting  the  calibrated  delays 
along  the  ground  station-to-spacecraft  line  of  sight. 

The  tropospheric  propagation  model  currently  implemented  in  the  geolocation  accuracy 
program  is  based  on  the  elevation  angle  dependent  model  developed  by  Hopfield  [1]. 

In  the  geolocation  program,  the  weight  function  is  used  as  the  inverse  covariance  of  noise 
Time  Difference  of  Arrival  (TDOA),  ephemeris  error,  altitude  error,  and  atmosphere 
error.  This  weight  function  is  used  to  evaluate  a  miss  distance  and  an  ellipse  uncertainty 
Semi-Major  Axis  and  Semi-Minor  Axis  (SMA  and  SMI).  The  current  geolocation  system 
used  the  Hopfield  atmosphere  propagation  delay  model  to  compute  the  atmosphere  error. 
At  below  0  degree  elevation  angle  from  the  horizon,  the  Hopfield  model  ignored  (cut  off) 
tropospheric  propagation  delay.  This  means  that  the  atmospheric  covariance  among  the 
weight  functions  is  neglected  in  computing  a  miss  distance  and  ellipse  containment  or  that 
the  measurements  for  TDOA  below  0  degree  from  the  horizon  are  considered  as  bad 
measurements.  Also,  the  Hopfield  model  uses  a  constant  refractivity  for  the  entire  world. 
In  this  report,  the  exponential  tropospheric  model,  which  provides  the  tropospheric 
propagation  delay,  even  at  negative  elevation  angles,  is  used  to  improve  the  geolocation 
accuracy  program.  The  exponential  tropospheric  model  is  derived  analytically  by  the 
mathematical  formulation  and  required  minimum  inputs  for  the  user.  This  exponential 
model  generates  the  refractivity  and  bending  error  adaptively  by  reading  real  time  surface 
meteorological  data  (temperature,  pressure,  dew  point  or  relative  humidity)  in  order  to 
process  more  accurate  results  of  tropospheric  effects  and  was  built  as  a  library  function  to 
ease  the  interface  with  any  other  applicable  system. 

In  this  report  the  chronology  of  the  tropospheric  software  and  database  development 
effort  is  presented.  Section  2  presents  the  improved  real-time  tropospheric  program  and 
the  binary  database.  Section  3  describes  the  data  extraction  and  processing  of  the 
procedure  for  FNL  data  used  to  generate  binary  format  databases  for  the  tropospheric 
program.  Section  4  shows  the  results  of  the  application  of  the  tropospheric  program  into 
the  geolocation  process.  The  geolocation  accuracy  of  the  Choi  model  improves 
approximately  50%  as  reflected  in  the  SMA  plots  and  miss  distance  calculations.  In 
Section  5,  conclusions  and  recommendations  derived  from  this  study  are  presented. 

Manuscript  approved  March  16, 2000. 
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2.  Improvement  on  Tropospheric  Model  and  Database 

The  first  version  of  tropospheric  propagation  model  (Choi  model)  code  [2]  is  designed 
as  a  subroutine  package  that  can  be  incorporated  into  an  application.  All  codes  are 
written  in  the  ANSI  standard  FORTRAN  77  with  no  machine  dependent  extensions 
except  the  VAX  system.  This  program  is  divided  into  the  Input/Output  (I/O)  data  file 
processing  and  the  main  calculation.  The  I/O  data  files  processing  program  reads  two 
data  files  in  ASCII  format.  The  first  contains  weather  data,  including  grid  number, 
latitude,  longitude,  surface  reffactivity,  pressure,  temperature,  and  undulation  of  the  geoid 
parameter.  The  second  contains  the  reference  height  coefficient  for  each  2.5°  x  2.5°  grid. 
This  part  of  the  program  invokes  a  database  for  the  specified  user  location  and  sends  it  to 
the  main  part  of  the  program  to  generate  the  tropospheric  propagation  effects.  The  main 
program  calculates  the  time  delay,  range  error,  and  angle  of  arrival  error  using  the  data 
files  from  the  I/O  part.  The  database  contains  the  ACSII  format  data  files  for  input  data 
and  output  data  files. 

The  new  version  of  the  tropospheric  program  has  been  converted  from  a  FORTRAN 
version  program  to  a  C  version  program.  This  new  C  version  program  was  upgraded  with 
a  new  algorithm,  and  was  built  as  a  library  function  to  easily  interface  with  the  other 
processors.  The  new  version  optimized  the  old  version  of  the  tropospheric  program  to 
adapt  to  the  real  time  operational  system,  and  converts  the  ASCII  database  format  to 
binary  format  to  reduce  processing  time  and  to  save  memory  size.  Figure  1  shows  a 
functional  block  diagram  for  the  new  tropospheric  propagation  program.  There  are  three 
major  sections  in  the  optimized  tropospheric  program:  1)  input  parameters  and 
determination  of  the  layer  height  and  reffactivity  profiles;  2)  the  apparent  angle 
computation  to  the  given  line  of  sight  (LOS)  elevation  angle;  and  3)  the  stratified  method 
to  calculate  the  time  delay,  range  error,  and  angle-of-arrival  error. 


Figure  1.  Functional  Block  Diagram  of  Exponential  Tropospheric  (Choi)  Model 
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In  this  new  version  of  the  tropospheric  model  program,  we  assume  that  the 
ionospheric  refractive  effects  are  negligible  above  400  MHz,  that  the  refractive  index 
gradient  exists  only  in  the  vertical  plane,  and  that  the  index  of  refraction  profile  can  be 
approximated  by  a  number  of  linear  segments  whose  thickness  is  small  compared  to  the 
earth’s  radius. 


2.1  Tropospheric  Propagation  Model  Program 

In  order  to  interface  with  any  real  time  processor,  the  tropospheric  delay  program 
must  be  implemented  as  a  library  routine  that  can  be  easily  maintained,  modified,  and 
ported  to  support  different  real  time  operational  platforms.  This  upgraded  C-version  of 
the  tropospheric  delay  program  was  developed  to  meet  real  time  data  processing 
requirements.  This  version  uses  a  macro  substitution  and  file  inclusion  as  a  pre¬ 
processor.  The  mathematical  functions  such  as  square  root,  cosine,  sine,  tangent,  and 
arctangent  for  the  degree  are  defined  as  macros.  The  variables  which  are  used  more  than 
twice  are  declared  as  “#defme”  statements.  The  input  database  sources  are  also  defined  as 
text  substitution  to  pass  the  library  function.  For  any  other  required  input  information 
and  functions  interfaced  with  library  functions,  “tropodelay.h”  is  used  with  “#include” 
statements.  In  order  to  reduce  function  generality  and  program  readability,  we  declare 
one  integer  array,  and  two  static  arrays  as  global  values.  All  functions  and  subroutines  in 
this  program  have  access  to  these  global  variables.  All  functions  in  the  library  should  be 
modularized.  This  means  that  a  program  should  be  broken  into  small  modules  with  one 
main  module  that  calls  the  others  through  the  proper  sequence.  The  “tropodelay”  function 
initially  connects  to  the  binary  database  to  access  the  user’s  desired  data.  This  function 
continuously  calculates  an  apparent  angle  to  the  line-of-sight  elevation  angle  from  the 
ground  station  to  a  target,  the  troposphere  delay  outputs  the  error  of  range,  the  angle  of 
arrival,  and  the  time  delay. 


2.1.1  Tropospheric  Height  and  Refraction  Profile 

The  index  of  refraction  in  the  troposphere  is  a  function  of  such  meteorological 
variables  as  temperature,  pressure,  and  water  vapor.  In  order  to  estimate  the  refractive 
effects  of  the  troposphere,  the  index  of  refraction  (n)  for  each  layer  height  must  be 
determined. 

The  new  C-version  program  assumes  that  the  troposphere  region  is  stratified  into  46 
different  height  and  constant  reffactivity  profiles  for  each  layer.  The  height  profiles  are 
increased  every  100  meters  between  0  and  1,000  meters,  every  500  meters  between  1,000 
and  10,000  meters,  and  every  1,000  meters  between  10,000  and  the  top  of  the 
troposphere  region  (27  km  from  the  Earth’s  surface).  The  main  reason  for  this  division  is 
based  on  the  fact  that  the  refractive  effects  are  small  above  10  km  from  the  surface  of  the 
Earth.  The  test  has  been  performed  in  order  to  validate  this  argument  of  46  layers  by 
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computing  range  and  angle  errors  on  different  environments  with  spatial,  temporal,  and 
geographical  variations  in  the  previous  study  for  more  than  1 00  sites  [2]. 

The  refractive  profiles  are  accomplished  by  a  combination  of  the  meteorological  data 
supplied  by  the  National  Oceanographic  and  Atmospheric  Administration  (NOAA)  and  the 
National  Weather  Service  (NWS)  and  the  exponential  tropospheric  propagation  model. 
The  exponential  tropospheric  propagation  model  (Choi  model)  that  the  NRL  has 
developed  is  used  to  estimate  the  refractivity  (N  =  106[n-l])  of  the  troposphere  region  (up 
to  27  Km  from  the  Sea  level)  based  on  the  height  of  46  different  profiles.  The  exponential 
tropospheric  propagation  profile  is  chosen  because  the  refractivity  of  the  observation  and 
prediction  is  very  close  to  real  weather  information  of  each  layer  in  the  troposphere 
region.  The  exponential  model  is  theoretically  based  on  the  surface  temperature,  pressure, 
and  relative  humidity.  Each  component  profile  starts  with  a  locally  observed  surface 
meteorological  value  and  decreases  to  the  effective  troposphere  height  (27  Km).  With  this 
assumption,  each  layer  refractivity  is  calculated  by  the  exponential  tropospheric  model 
program.  Empirical  weather  data  is  used  in  the  calculation  of  refractivity  for  each  of  the 
layer  profiles  when  local  weather  information  is  not  available. 

Figures  2  and  3  show  the  different  time  delays  for  different  azimuth  angles  for 
February  and  August  between  the  original  ECM  14  layer  data  and  the  model  defined  by  46 
layer,  respectively,  for  0  degree  elevation  angle  from  the  ground  to  a  target  for  8  different 
locations  (refer  to  Appendix).  In  these  figures,  the  different  time  delays  for  azimuth  angle 
(0  ~  360  degree)  are  within  ±  0.05  nanosecond  for  8  different  sites.  Figures  4  and  5 
present  time  delay  differences  between  original  ECM  14  layer  data  and  model  defined  by 
29  layer.  Figures  6  and  7  show  the  time  delay  difference  between  ECM  14  layer  data  and 
the  model  defined  by  14  layer.  From  these  figures,  the  conclusion  has  been  made  to  use 
46  different  layer  profiles  and  heights  because  of  their  accuracy  and  closeness  to  the 
analytical  result. 


2.1.2  Apparent  Angle  Computation 

Since  the  refractive  index  is  not  the  same  everywhere,  the  trajectory  tends  to  be  curved 
consequently.  There  is  a  difference  between  the  real  direction  of  the  source  and  the 
apparent  direction  from  which  the  propagation  waves  arrive.  We  have  already  seen,  in  the 
case  of  the  troposphere,  that  in  a  spherically  symmetrical  medium  the  radius  of  curvature 
of  the  trajectory  is  proportional  to  the  vertical  gradient  of  the  refractive  index.  The  first 
version  of  the  tropospheric  program  was  used  as  an  iterative  approach  to  compute  an 
apparent  angle  given  the  LOS  (line-of-sight)  elevation  angle  from  the  ground  to  satellite. 
Figure  8  presents  the  apparent  angle  computation  algorithm  using  an  iterative  approach. 
If  the  angle  error  difference  between  a  LOS  angle  and  an  estimated  apparent  angle  is 
within  a  given  threshold,  the  convergence  is  reached  and  the  true  incoming  elevation  angle 
is  now  known.  An  initial  estimate  of  the  elevation  angle  for  the  actual  signal  is  obtained 
by  adding  a  different  angle  of  arrival  error  (AOA)  to  the  LOS  elevation  angle 
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Figure  2.  Difference  of  Time  Delay  (ECM  14  layer  and  model  46  layer)  for  February 


Figure  3.  Difference  of  Time  Delay  (ECM  14  layer  and  model  46  layer)  for  August 
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Figure  4.  Difference  of  Time  Delay  (ECM  14  layer  and  model  29  layer)  for  February 


Figure  5.  Difference  Time  Delay  (ECM  14  layer  and  model  29  layer)  for  August 
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Figure  6.  Difference  of  Time  Delay  (ECM  14  layer  and  model  14  layer)  for  February 
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Figure  7.  Difference  of  Time  Delay  (  ECM  14  layer  and  model  14  layer)  for  August 


Figure  8.  Iteration  Method  for  Apparent  Angle 

Snell’s  law  for  spherical  geometry  and  the  stratified  layer  method  with  the  refractivity 
height  profile  are  used  to  compute  the  signal  path  at  each  location  in  the  troposphere.  If 
the  deviation  estimated  apparent  angle  from  the  true  LOS  elevation  angle  is  above  a 
certain  threshold,  the  process  is  repeated  until  convergence  is  reached.  Above  10  degrees 
LOS  elevation  angle,  the  iteration  loop  is  not  necessary  to  reach  the  convergence 
condition.  Below  the  10  degree  LOS  elevation  angle,  more  than  one  iteration  loop  is 
needed  to  satisfy  the  convergence  condition.  Below  a  10  degree  LOS  elevation  angle,  at 
least  2  iteration  loops  are  specifically  requested.  An  iterative  approach  to  compute  the 
apparent  angle  for  the  refractive  effect  is  not  adaptable  to  the  real-time  processor  because 
the  processing  time  of  the  system  is  too  long. 

In  order  to  reduce  processing  time  of  the  real  time  processor,  we  have  analyzed  and 
developed  a  new  algorithm  for  the  apparent  angle  using  the  given  input  parameters  of  the 
troposphere  propagation  delay  model.  Figure  9  shows  the  angle  error  versus  refractivity 
for  all  seasons  to  different  LOS  elevation  angles.  In  this  figure,  the  apparent  angle  can  be 
predicted  with  the  LOS  elevation  angle  and  surface  refractivity,  which  is  observed  from 
simple  surface  measurements  of  the  common  meteorological  elements  of  temperature, 
pressure,  and  relative  humidity.  Each  linear  equation  is  associated  with  the  LOS  elevation 
angle  and  the  surface  refractivity.  The  new  statistical  apparent  angle  computation  is  based 
on  the  slope  and  the  intercept  component  of  each  linear  equation.  The  apparent 
polynomial  equation,  which  is  combined  with  the  slope  and  intercept  component,  is  a 
function  of  the  surface  refractivity  and  the  LOS  elevation  angle.  This  new  algorithm  was 
inserted  in  the  new  improved  tropospheric  propagation  model  program.  Figure  10  is 
shown  for  the  difference  of  time  delay  between  the  iterative  approach  and  the  statistical 
approach  for  the  verification  of  the  new  algorithm  for  apparent  angle  computation  of  8 
different  locations.  In  these  figures,  the  absolute  value  of  angle  differences  are  within  0.3 
nanosecond. 
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Figure  10.  Time  Delay  Difference  between  Iteration  and  Statistical  Approach 


2.2  Procedure  to  Reduce  Database  and  Extract  Binary  Data 

The  previous  C-version  of  the  tropospheric  delay  program  only  supported  ASCII  data 
files,  which  read  sequentially  within  the  inner  loop  when  data  is  requested  from  a  file.  This 
process  adds  an  additional  overhead  to  the  program  execution  time.  This  ASCII  format 
database  needs  a  huge  storage  space  for  high-resolution  global  data  such  as  the  MRF 
hourly  and  diurnal  database.  We  found  that  the  binary  data  format  could  be  processed 
faster  than  ASCII  format  data.  When  we  switched  raw  data  to  the  binary  format,  we  used 
a  simple  compression  technique  that  converted  floating  point  values  (4  bytes)  to  integer 
values  (2  bytes).  The  previous  data  access  routine  was  invoked  by  the  data  based  on  the 
grid  number  from  the  database.  In  this  new  version,  we  used  the  record  number  instead  of 
the  grid  number  for  direct  access  to  the  data  and  binary  tree  linked  list  approach.  When 
the  user  selects  input  data  information,  the  new  program  calculates  the  offset  value  for  the 
record  instead  of  going  back  to  the  grid  number  routine.  Therefore,  the  new  binary  data 
format  does  not  need  to  keep  the  grid  number.  This  process  saves  42  Kbytes  for  the  ECM 
database,  168  Kbytes  for  HIRAS,  and  61  Mbytes  for  MRF.  Table  1  shows  a  comparison 
of  the  database  size  between  ASCII  and  the  binary  format.  The  size  of  the  database  is 
reduced  by  up  to  approximately  50%  for  each  database  through  the  new  program.  The 
modified  version  will  read  binary  files,  and  directly  access  to  the  data.  The  inner  loop  was 
then  removed  from  file  readings. 
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Table  1.  Comparison  of  ASCII  and  Binary  Format  Size 


ASCII  Format 

Binary  Format 

ECM 

2  Mbytes 

1  Mbytes 

HIRAS 

5.4  Mbytes 

2.5  Mbytes 

MRF 
(7  Years) 

1.3  Gigabytes 

636  Mbytes 

FNL 

(Monthly) 

384  Mbytes 

109  Mbytes 

FNL 

(Yearly) 

4.7  Gigabytes 

1.15  Gigabytes 

ECM  Binary  Format  -  ‘Yearly  data’  contains  10,368  records  with  72  bytes  for 
each  record  as  shown  in  Figure  11.  The  first  24  bytes  of  the  record  contain  surface 
reff activity  data  for  12  months.  The  last  48  bytes  of  the  record  contain  the  reference 
height  data  for  the  12  months.  Each  monthly  surface  refractivity  data  must  be  read  as  2 
bytes  unsigned  integer  or  short.  The  first  2  bytes  of  the  24  bytes  refer  to  the  first  month, 
the  next  2  bytes  of  the  24  bytes  refer  to  the  second  month,  and  so  forth.  To  retrieve  the 
actual  refractivity  value,  you  must  extract  the  first  2  bytes  of  data  related  to  the  offset 
value,  and  then  convert  the  two-byte  integer  value  to  an  actual  refractivity  floating  point 
divided  by  100.0.  Each  monthly  reference  height  data  must  be  read  as  4  bytes  long.  The 
first  4  bytes  of  the  48  bytes  refer  to  the  first  month,  the  next  4  bytes  of  the  48  bytes  refer 
to  the  second  month,  and  so  forth.  To  retrieve  the  actual  reference  height  value,  you  must 
convert  the  extracted  value  to  a  floating  point  value  and  divide  it  by  100.0. 

HIRAS  Binary  Format  -  ‘Monthly  6  hourly  data’.  The  HIRAS  binary  file  is 
divided  into  monthly  6-hour  blocks  of  data.  Each  monthly  block  of  data  consists  of 
124,416  bytes,  and  contains  a  6-hour  block  of  data  for  all  10,368  2.5°  x  2.5°  grids.  Each 
grid  block  of  data  is  12  bytes  in  size  as  shown  in  Figure  12.  The  first  8  bytes  contain  6- 
hour  surface  refractivity  data,  and  the  last  4  bytes  contain  the  reference  height  data.  The 
0000-hour  would  be  the  first  2  bytes,  the  0600  hour  would  be  the  next  2  bytes,  and  so 
forth.  Surface  refractivity  data  is  read  as  2  bytes  unsigned  integer  or  short.  The  extracting 
method  of  the  actual  value  is  the  same  as  in  the  case  of  the  ECM  database. 


11 


Figure  11.  ECM  Binary  Format  Structure 


Figure  12.  HIRAS  Binary  Format  Sructure 
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MRF  Binary  Format  -  ‘Daily  6-hour  data’.  The  MRF  binary  file  is  divided  into 
yearly  binary  files.  Each  yearly  file  contains  1 2  months  of  6-hour  daily  data.  Each  6-hour 
daily  block  of  data  is  248,832  bytes  in  size,  and  contains  all  10,368  2.5°  x  2.5°  grids  for 
that  particular  day.  Each  grid  block  of  data  is  24  bytes  in  size  as  shown  in  Figure  13.  The 
first  8  bytes  contain  6-hour  surface  refractivity  data,  and  the  last  16  bytes  contain  the  6- 
hour  reference  height  data.  The  0000-hour  would  be  the  first  2  bytes,  the  0600  hour 
would  be  the  next  2  bytes,  and  so  forth.  Surface  refractivity  data  is  read  as  2  bytes 
unsigned  integer  or  short.  The  last  16  bytes  contain  the  6-hour  as  it  is  reference  height 
data.  The  0000  hour  would  be  the  first  4  bytes  out  of  the  16  bytes,  the  0600  hour  would 
be  the  next  4  bytes  out  of  the  16  bytes,  and  so  forth.  Reference  height  data  must  be  read 
as  4  bytes  long.  The  actual  value  is  extracted  in  the  same  method  as  for  the  case  of  ECM. 


Figure  13.  MRF  Binary  Format  Structure 
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2.3  Performance  Analysis  and  Simulation  Results 

The  improved  tropospheric  exponential  model  with  binary  database  has  been 
examined  to  show  the  advantages  and  improvements  of  the  new  approach.  Figures  14, 
15,  and  16  show  global  time  delay  errors  for  1  degree,  0  degree,  and  -1  degree  elevation 
angle  from  the  horizon  during  August  based  on  a  10-year  average  of  climatology  data 
from  the  ECMWF.  Time  delays  vary  from  180  nanoseconds  in  the  shorter  range  to  620 
nanoseconds  in  the  longer  range.  This  error  reflects  a  180  ft  range  error  for  the  ideal 
cases  without  consideration  of  other  atmospheric  effects.  Figures  17,  18,  and  19  show 
contour  plots  of  global  time  delay  errors  for  the  same  elevation  angle  with  the  FNL  data 
during  August  15,  1997  at  00:00  Zulu  time  based  on  diurnal  meteorological  data  from 
NOAA.  Those  figures  show  the  different  time  delays  between  average  database  and  real 
time  database  for  the  entire  world.  Table  2  shows  the  time  delay  comparison  between  a 
10-year  average  empirical  database  (ECM)  and  real  time  data  for  11  different  locations 
from  Figures  14  through  19.  This  table  provides  the  time  delay  difference  from  the 
negative  3-degree  elevation  angle  to  the  10-degree  elevation  angle.  The  variation  of  time 
delay  is  noticeable.  For  below  3-degree  elevation  angle  from  the  horizon,  the  time  delay 
difference  dramatically  increases  in  the  order  of  magnitude.  The  result  of  Table  2  indicates 
that  the  choice  of  the  real-time  data  dominates  the  performance  of  the  atmospheric 
propagation  delay.  This  table  shows  that  the  time  delays  above  a  5-degree  elevation  angle 
are  not  much  different  from  those  below  a  3-degree  elevation  angle.  Table  3  shows  the 
angle  of  arrival  error  comparison  from  10  degrees  to  negative  3  degrees  between  the 
exponential  model  and  theHopfield  model  using  the  real  time  FNL  data  August  1,  1999 
data  for  8  different  locations.  In  this  table,  we  know  that  the  exponential  model 
compensated  more  than  50%  for  any  elevation  angle.  This  wide  error  causes  significant 
problems  in  satellite  communication  link  analysis  and  antenna  operator.  The  time  delay 
and  angle-of-arrival  error  produce  a  geolocation  error  of  tens  of  nautical  miles.  If  these 
errors  are  calibrated  or  corrected  on  time  for  the  specific  mission,  valuable  time  and  cost 
can  be  saved,  and  the  accuracy  will  be  dramatically  increased. 

The  improved  tropospheric  exponential  model  with  binary  database  reduces  the 
processing  time  when  we  calculated  time  delay,  range  error,  and  angle  of  arrival  error. 
Table  4  shows  the  average  processing  time  for  different  layers  with  and  without  grid 
crossing.  The  processing  times  were  computed  by  executing  an  iteration  loop  3,600 
times.  The  number  in  parenthesis  represents  the  processing  time  of  one  iteration  of  the 
tropospheric  program.  The  processing  time  of  binary  format  is  faster  than  ASCII  format 
in  this  table. 


3.  Data  Extraction  and  Processing  Procedure 

From  March  1999  to  date,  we  have  processed  the  FNL  and  MRF  database  (1997  - 
present).  We  are  currently  developing  an  automatic  and  user  friendly  data  processing 
software  tool  to  produce  the  database  for  the  new  tropospheric  program.  The  FNL  data 


14 


90 


. 

•:: -v.s •  > -': '•  .-^si-ti c«ss-~ .-'  '*••; 5*Sf ,•> ^-£L?lr?J?#&Sv: 


•:-n^: 


tSfcciii 

"00 


,^iiiii 

MltilSlMIS! 


«»^fe^qaj^fS8s 


2 

E 

a 

S 

2 

ECM  Time  Delay  Contour  Map  for  August  .  Truezenang  =  91.0  deg 


Figure  16.  ECM  Time  Delay  Contour  Map  for  August  (elevation  angle  =  -1.0  deg) 


=|\| L Time  Delay  Contour  Map  for  August  15,1  997  (OOOO  hour)  .  T ruezenang  =69  0  deg 


Figure  17.  FNL  Time  Delay  Contour  Map  for  August  15,  1997  (elevation  =1.0  deg) 
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Table  2.  Time  Delay  Difference  between  ECM  and  FNL  for  August  1997 


Table  3.  Angle  Error  Comparison  August  1,  1999  for  8  sites 
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file  contains  hemispheric  129  by  129  stereographic  grids  and  data  on  one  hemisphere  for 
approximately  half  a  month;  days  one  through  15,  and  16  through  the  end  of  month. 
The  raw  FNL  data  are  stored  in  the  unpacked  binary  form  with  the  header  record.  Each 
data  record  was  composed  of  a  50-character  header  in  ASCII  format,  followed  by  the 
binary  packed  data.  To  read  the  header  record  and  extract  binary  data,  one  should 
decode  the  binary  index  record.  To  generate  the  FNL  database  of  1  degree  by  1  degree 
high-resolution  world  data  for  the  tropospheric  program,  one  should  reprocess  raw  FNL 
data.  Figure  20  shows  block  diagram  of  the  on-line  data  processing  by  connecting  to  the 
NOAA  database  to  generate  the  database  for  a  tropospheric  exponential  model  program. 


Figure  20  On-line  FNL  Data  Processing  Diagram 


3.1  Unpacking  FNL  Database 

The  FNL  data  processing  begins  with  the  unpacking  of  the  field  and  printing  out 
the  grid  point  which  corresponds  with  the  latitude  and  longitude.  The  program 
"unpack  fill"  performs  the  unpacking  and  generates  the  surface  refractivity  data  from  the 
extracted  data.  If  the  whole  data  set  is  missing,  the  program  creates  a  "Missing.dat"  file. 
This  file  contains  information  on  what  data  is  missing,  and  returns  the  number  of  bytes  to 
be  read.  In  this  routine,  the  input  data  file  is  the  raw  FNL  data  file  from  NOAA.  This 
routine  passes  the  record  number  and  the  temporary  output  file  name  to  store  unpacked 
data  for  each  record  number  as  input  parameters.  The  array  filename  is  used  to  generate 
the  name  of  the  selected  record.  The  following  is  an  example  of  a  relation  between  the 
record  number  and  the  filename. 
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rec  1  =  hour  0000  day  1  -  fnl  month  O  l_year  00. final 

rec  2  =  hour  0600  day  1  -  fhl_month_01_year_06.  final 

rec  3  =  hour  1200  day  1 
rec  4  =  hour  1 800  day  1 

rec  5  =  hour  0000  day  2  -  fhl_month_02_year_00.  final 

rec  6  =  hour  0600  day  2  -  fnl_month_02_year_06.final 


According  to  the  record  number,  the  routine  saves  data  in  the  memory  buffer.  Then, 
reads  the  raw  unpacked  data  file  to  create  the  output  data  file.  One  buffer  contains  14 
different  levels  of  data  based  on  different  pressure  from  surface  to  20  mb  at 
approximately  27  km  height.  Using  the  iteration  loop,  each  layer  data  is  extracted  from 
the  buffer.  The  output  file  produces  a  3,261,636  byte  binary  file  that  contains  a  129  by 
129  stereographic  grid  of  the  unpacked  atmospheric  data.  Each  grid  is  14  bytes  in  size  as 
shown  in  Figure  21.  The  first  4  bytes  contain  pressure  data.  The  next  2  bytes  contain 
temperature  data.  The  next  2  bytes  contain  relative  humidity  data.  The  next  2  bytes 
contain  surface  refractivity  data.  The  last  4  bytes  contain  reference  height  data. 


Figure  21.  Unpacking  Processing  Diagram 
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3.2  Merging  Processing 

There  is  no  latitude,  longitude,  and  surface  height  information  in  the  new  unpacked 
FNL  data  file.  We  connect  to  the  NOAA  Air  Resource  Laboratory  in  Maryland  to  get 
the  location  information  for  129-by-129  polar  stereographic  grids.  They  send  two 
northern  and  southern  hemispheric  data  sets.  We  generate  a  new  file  that  merges  the 
129-by-129  latitude  and  longitude  data  with  the  atmospheric  data.  This  merging  program 
creates  new  temporary  binary  data  files  that  merge  latitude  and  longitude  for  each  grid. 
Each  grid  block  size  of  a  new  binary  data  file  is  22  bytes  as  shown  in  Figure  22.  The 
first  4  bytes  are  for  latitude  data  and  the  next  4  bytes  for  longitude  data.  The  next  4  bytes 
are  for  pressure  data  and  the  next  2  bytes  for  temperature  data.  The  next  2  bytes  are  for 
relative  humidity  data  and  the  next  2  bytes  for  surface  reffactivity  data.  The  last  4  bytes 
contain  reference  height  data. 


Figure.  22  Merge  Data  Processing  Diagram 
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3.3  Calculating  Reference  Height 


The  tropospheric  exponential  model  program  needs  a  reference  height  for  each 
grid  as  a  scale  factor.  The  newly  created  file  contains  latitude,  longitude,  and 
atmospheric  data  to  generate  a  new  binary  file  with  reference  height  derived  from 
atmospheric  data.  To  determine  the  reference  height  for  each  grid,  a  simple  linear 
interpolation  method  has  been  adopted.  Given  the  refractivity/height  pairs  h,/n,  and 
h2/n2,  we  desire  to  determine  the  reference  height  h0  at  reffactivity  n0  (see  Figure  23). 


Note  that, 


And 


log  h2  =  mlog«2  +  b 


log/z,  =  mlogH,  +b 


Solving  for  m  and  b, 


m-  MV'*,] 

l°g  [«,/«2] 


b  =  logk,  -  log /z, 


(1) 

(2) 


(3) 

(4) 


The  reference  height  h^  at  the  reffactivity  r^  is  then  found  from 

log/zo  =  /wlog«0  +b  (5) 


Figure  23.  Interpolation  Method 
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The  return  value  of  this  routine  points  to  the  new  output  file.  Each  grid  block  size 
of  the  output  file  is  22  bytes.  The  first  4  bytes  are  for  latitude  data  and  the  next  4  bytes 
for  longitude  data.  The  next  4  bytes  are  for  pressure  data  and  the  next  2  bytes  for 
temperature  data.  The  next  2  bytes  are  for  relative  humidity  data  and  the  next  2  bytes  for 
surface  refractivity  data.  The  last  4  bytes  contain  reference  height  data.  Figure  24  shows 
the  diagram  of  calculating  reference  height. 


Figure  24.  Calculating  Reference  Height  Diagram 
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3.4  Averaging  Atmospheric  Data 

There  are  more  than  2  information  blocks  at  the  same  grid  when  the  FNL  data  is 
reprocessed  to  1  degree  by  1  degree  for  the  tropospheric  program.  We  use  the  newly 
created  and  averaged  atmospheric  data  when  grid  numbers  of  the  same  value  intersect.  A 
new  binary  file  is  created  which  contains  the  calculated  grid  numbers  and  averaged 
atmospheric  data.  All  data  in  the  file  is  sorted  in  ascending  order  according  to  the  grid 
number.  The  binary  file  data  is  sorted  by  using  a  binary  tree  sorting  algorithm.  A  binary 
tree  is  a  hierarchical  structure  based  on  a  modified  form  of  a  linked  list,  which  is  request- 
and-release  memory  ffom/to  the  operating  system  as  required  by  the  program  at  a  run¬ 
time.  Each  node  in  the  binary  tree  maintains  two  further  node  pointers,  known  as  ‘left’ 
and  ‘right’,  which  point  to  the  left  and  the  right  sub-trees  of  that  node.  Absent  sub-trees 
are  marked  by  a  NULL  valued  left  or  right  within  their  corresponding  node.  A  node  with 
neither  left  nor  right  sub-tree  is  known  as  a  leaf  node,  while  the  node  at  the  top  of  the 
tree  is  called  the  root  node.  This  binary  tree  sorting  algorithm  has  the  ability  to  rapidly 
search  and  sort  the  data.  Also,  it  is  easy  to  insert  and  delete  data  into  memory.  The  data 
structure  of  this  binary  tree  is 

Structure  BinNode  { 

Int  item; 

Struct  BinNode  *left; 

Struct  BinNode  *right;  } 


The  integer  item  presents  the  grid  number.  The  grid  number  is  greater  than  the  grid 
number  of  the  other  node,  then  it  goes  to  the  right  binary  node.  It  does  not  go  to  the  left 
node  to  sort  the  grid  number  for  random  data  file.  According  to  the  sorted  grid  number, 
the  data  in  the  memory  buffer  copies  into  the  record  array  to  create  an  output  file.  The 
output  file  of  this  routine  contains  grid  number  data  instead  of  latitude  and  longitude  data 
for  each  grid.  Each  grid  block  size  of  output  file  is  16  bytes.  The  first  2  bytes  are  for 
grid  number  data  and  the  next  4  bytes  for  pressure  data.  The  next  2  bytes  are  for 
temperature  data  and  the  next  2  bytes  for  relative  humidity  data.  The  next  2  bytes  are  for 
surface  reffactivity  data  and  the  last  4  bytes  for  reference  height  data.  Figure  25  shows 
the  average  stratified  method  data  processing  diagram. 


3.5  Missing  Data 

Once  data  has  gone  through  the  averaging  atmospheric  data  process,  it  is  then  highly 
probable  that  there  may  be  some  missing  data.  All  64,800  grid  numbers  may  not  have 
been  generated.  To  fill  these  missing  grid  numbers  with  their  corresponding  linear 
interpolated  atmospheric  data,  linear  interpolation  techniques  were  applied.  This 
interpolation  routine  checks  and  generates  a  final  binary  file  by  the  grid  number  with 
their  corresponding  linear  interpolated  atmospheric  data.  In  this  routine,  we  used  a 
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Figure  25.  Average  Stratified  Method  Data  Processing  Diagram 
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double  linked-list  algorithm  for  the  missing  data  point.  The  first  two  nodes  (head  and 
tail)  are  initialized  and  invoked  into  node.  If  the  difference  of  the  grid  number  into  the 
head  and  tail  node  is  greater  than  one,  missing  data  points  exist  in  two  point  nodes.  This 
routine  creates  nodes  as  much  as  the  difference  between  the  head  and  tail  grid  number 
and  fill  missing  data  using  the  linear  interpolation  method  into  the  newly  created  node. 
The  final  grid  block  size  of  each  output  file  is  16  bytes.  The  first  2  bytes  are  for  the  grid 
number  data  and  the  next  4  bytes  for  pressure  data.  The  next  2  bytes  are  for  temperature 
data  and  the  next  2  bytes  for  relative  humidity  data.  The  next  2  bytes  are  for  surface 
reff activity  data  and  the  last  4  bytes  for  reference  height  data.  Figure  26  shows  the 
additional  missing  data  processing  diagram. 
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Figure  26.  Additional  Missing  Data  Processing  Diagram 
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3.6  Automatic  Process 


A  shell  script  was  developed  for  the  execution  of  each  routine  to  process  raw  NOAA 
FNL  binary  files  automatically.  There  are  two  input  options  to  execute  the  shell  script. 
One  automatically  calculates  how  many  records  are  available  to  generate  the  final  FNL 
data  file.  The  other  processes  the  particular  record  only.  This  routine  processes  the 
routine  3.1  through  routine  3.5  automatically  with  the  input  record  parameter  to  create 
the  final  output  binary  data  files. 


4.  Geolocation  Accuracy  Measure  with  Tropospheric  Delay 

The  geolocation  of  a  radio  transmitter  has  a  wide  variety  of  applications.  One 
example  is  to  locate  an  emergency  radio  source  for  ship  or  aircraft  rescue.  Another  is  to 
find  the  position  of  an  interference  source  that  interferes  with  satellite  operations. 

The  geolocation  measure  of  the  radio  transmitter  is  based  mainly  on  techniques 
which  relate  to  one  or  a  combination  of  frequency  and  time  information.  The  set  of 
delays,  or  time  difference  of  arrival  (TDOA),  associated  with  each  of  the  source  wave 
fronts  between  the  receivers  can  be  used  to  determine  the  source  location  by  solving  a 
nonlinear  or  linear  system  equation.  In  this  report,  we  use  TDOA  measurements  in  which 
a  propagation  time  for  the  uplink  signal  through  a  particular  satellite  is  compared  with  an 
adjacent  satellite.  Given  that  the  position  of  two  spacecraft  relative  to  the  receiving 
station  is  precisely  known,  the  difference  between  the  time  of  arrival  over  the  two  different 
paths  will  localize  the  uplink  transmitter  location  to  a  curve  on  the  surface  of  the  Earth. 
Since  the  TDOA  technique  described  only  localizes  the  position  of  the  uplink  transmitter 
to  a  curve  on  the  surface  of  the  Earth,  some  additional  techniques  are  required  to  resolve 
the  ambiguity.  Initially,  it  was  thought  that  the  ambiguity  could  be  resolved  by  making  the 
third  spacecraft  to  form  another  signal  path  pair.  Then,  the  position  of  the  uplink  in 
question  could  be  determined  by  finding  the  intersections  of  the  two  resulting  arcs.  Figure 
27  presents  a  two  dimensional  view  of  geometry  pertinent  to  the  TDOA  technique.  The 
position  of  the  two  satellites,  the  location  of  the  TDOA  site,  and  the  distance  from  the 
satellites  to  the  TDOA  receiver  site  are  known.  The  unknowns  are  the  two  remaining 
distances  from  each  of  the  satellites  to  the  unknown  uplink  station. 

If  the  receiver  site  has  two  antennas  connected  to  two  receiving  systems  which  can 
simultaneously  monitor  the  signal  through  each  satellite,  a  receive  time  delay  between  the 
two  paths  can  be  measured.  Given  this  delay  along  with  the  knowledge  of  the  speed  of 
the  signal  propagation  through  space  and  the  equivalent  distance  of  any  unnoticeable 
signal  processing  time,  the  difference  in  path  distances  may  be  deduced.  Subtracting  out 
the  known  downlink  distance  and  any  other  system  offsets  from  the  satellites  to  the  TDOA 
receiver  site  leaves  the  remaining  differential  distance  between  the  two  unknown  uplink 
distances. 
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Figure  27.  Two-Dimensional  TDOA  Geometry  (Isochrone  Intersection) 


4.1  Covariance  and  SMA/SMI  Estimation 

The  components  of  an  n-dimensional  vector  x  that  is  to  be  estimated  are  the 
position. 


x  =  [x,,  x2,  x„]  (6) 

A  set  of  N  measurements  yj,  i=  1,  2,  ...  n  is  collected  at  various  positions.  In  the  absence 
of  random  measurement  errors,  y;  is  equal  to  a  known  function  f(x).  In  the  presence  of 
additive  errors, 


yi  =  fi(x)  +  Vi,  i  =  1,2,  ...,  n  (7) 

These  n-equations  can  be  written  as  a  single  equation  for  n-dimensional  column  vectors; 

y  =  f(x)  +  v  (8) 
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where  v  is  measurement  noise,  the  estimate  of  x  is  a  matter  of  minimizing  a  cost  function 
such  as  the  weighted  sum  of  squared  residuals, 

Q(x)  =  e(x)‘  W  e(x)  (9) 

where  e(x)  =  y  -  f(x)  and  t  is  the  transpose  of  the  vector  e(x). 

The  weight  matrix  is  the  inverse  of  the  covariance  matrix  of  noise  measurements.  In 
general,  f(x)  is  a  nonlinear  vector  function.  To  determine  a  reasonably  simple  estimator, 
f(x)  can  be  linearized  by  expanding  it  in  a  Taylor  series  about  a  reference  point  specified 
by  the  vector  x0  and  retaining  the  first  two  terms;  that  is,  we  use 

f(x)  =  f(x0)  +  A  (x  -  xo  )  (10) 

Where  x  and  x0  are  n  x  1  column  vectors  and  A  is  the  n  x  n  matrix  of  derivatives  evaluated 
at  x0: 
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Each  row  of  this  matrix  is  the  gradient  vector  one  of  the  components  of  function  f(x).  The 
vector  xo  could  be  an  estimate  of  x  determined  from  a  previous  iteration  of  the  estimation 
procedure  or  based  on  prior  information.  It  is  assumed  in  the  subsequent  analysis  that  xo 
is  sufficiently  close  to  x  that  the  equation  f(x)  =  f(x<j)  +  A(x-Xo)  is  an  accurate 
approximation.  Combining  equation  (9)  and  equation  (10)  gives 

Q(x)  =  (e,  -  Axf'W^ei  -  Ax)  (12) 


where  ei  =  e  -  f(xo)  +  Axo- 

To  determine  the  necessary  condition  for  the  estimator  x  that  minimizes  Q(x),  we  calculate 
the  gradient  of  Q(x),  defined  by 


dQ 

dxl 


dQ 


3x2 


dQ 


(13) 
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And  then  solve  for  x  such  that 


^-Q(*)  =  °  (14) 


From  its  definition,  W  is  a  symmetric  matrix.  Therefore 

d  A 

—Q(x)=2A'WAx-2A'W  (15) 

We  assume  that  the  matrix  A‘WA  is  non-singular.  Thus  the  solution  is 


x  -  (A'WAy'A'Wei  =  x0+(A'WA)~' A' We  (16) 

For  all  cases,  the  covariance  matrix  [A‘WA]'’  at  the  final  step  gives  the  statistics  of  final 
position  estimates.  For  the  estimator  of  a  two-dimensional  vector,  such  as  position 
coordinates  on  the  surface  of  the  Earth,  the  bivariate  covariance  matrix  can  be  expressed 
as 


cr, 


P  = 


a 


12 


(17) 


and  with  the  assumption  that  the  error  distribution  is  a  Gaussian,  the  parameter  of  error 
ellipse  with  a  =  semi-major  axis  and  b  =  semi-minor  axis  is  determined  by 
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4.2  Geolocation  Measure  using  TDOA 

The  geolocation  program  uses  time  difference  of  arrival  (TDOA)  measurements  as 
the  measurement  parameter  for  its  estimation  process.  Thus,  the  estimation  problem 
becomes  a  problem  of  determining  the  functional  relationship  between  TDOA  and  the 
estimation  parameters.  Given  the  range  distance  from  the  emitter  to  the  spacecraft,  R(t) 
and  the  atmospheric  delay  datm(t),  the  TDOA  is  given  by 


x  (x,t)  =  (1/c)  R(t)  +  datm(t)  +  e 


(21) 


R(t)  =  ||  r(t)  ||  =  (r*  r) '/:  (22) 

r  =  s  -  1 

with  s  representing  the  position  of  spacecraft,  and  1  the  emitter  location.  Then,  the 
calculated  TDOA  between  spacecraft  1  and  2  is  given  by 

fl20O  =  X2(t2)  -  Tl(tl) 

=  (l/C)[R2(t2)  -  Rl(tl)]  +  d  atm,2(t2)  -  d  atm,l(tl)  (23) 


4.3  Analysis  Result 


The  TDOA  analysis  for  an  emitter  geolocation  is  required  to  validate  the  accuracy 
of  the  model  due  to  the  change  of  velocity  for  the  radio  frequency  (RF)  wave  with  respect 
to  the  refractivity  variation  in  the  troposphere  region.  The  precision  geolocation 
algorithms  of  the  current  Ramrod  system  use  an  approximation  of  the  Hopfield  model  that 
produces  significant  errors  in  computing  the  time  delay  at  low  elevation  angles  (less  than  3 
degrees)  The  Hopfield  model  applied  zero  time  delay  at  negative  elevation  angle  to  the 
spacecraft.  The  exponential  model  (Choi  model)  has  used  the  historical  climatology 
weather  data  and  an  improved  algorithm  to  provide  a  more  accurate  time  delay  profile. 
The  time  delay  profile  of  the  exponential  model  (Choi  model)  would  provide  better 
location  and  reduce  the  uncertainty  associated  with  those  inferred  emitter  positions.  The 
geolocation  program  has  been  developed  to  allow  the  operator  to  choose  which  of  the  two 
troposphere  models  (Hopfield  or  Choi  model)  should  be  adopted  for  the  process  of  a 
particular  signal  set.  A  100  nanoseconds  difference  between  the  Hopfield  and  Choi  model 
at  low  elevation  angles  suggests  that  a  significant  correction  may  be  able  to  be  made  in  the 
low  elevation  TDOA.  The  Choi  model  used  temporal  and  meteorological  data.  We  used 
different  atmospheric  factor,  covariance  time  delay,  and  atmospheric  delay  instead  of 
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constant  values  as  for  the  Hopfield  model  when  applied  to  the  current  Ramrod  system. 
We  observed  that  the  Choi  model  improved  the  miss  distance  and  the  uncertainty 
approximately  50%  in  comparison  with  the  current  Hopfield  model.  This  report  covered 
the  results  of  the  emitter  geolocation  data  sets  for  six  different  events  collected  in  January 
21,  1998.  Figure  28  depicts  the  collected  benchmark  data  format. 


Benchmark 

#1 

Benchmark 

#n 

Benchmark 

tn 

Benchmark 

Total  2034  benchmarks 

Figure  28.  Collected  Benchmark  Data  Format 


These  events  provide  a  typical  collection  of  signals  on  benchmark  channels.  The 
data  sets  associated  with  each  event  contain  the  scan  files  and  signal  files  associated  with 
the  events.  The  scan  data  included  the  initial  location  of  geolocation  program,  TOA  of 
each  spacecraft,  and  TDOA  associated  with  each  pair  of  the  spacecraft.  The  data  set  is 
applied  to  the  geolocation  program  for  the  results  of  miss  distance  and  the  uncertainty  of 
the  benchmark  emitter  position.  The  actual  benchmark  emitter  location  is  used  as  an 
initial  location  instead  of  the  estimated  location.  The  output  (Miss  distance,  SMA,  and 
SMI)  of  the  geolocation  program  is  iterated  only  once  for  the  above  scan  data  set.  The 
geolocation  program  read  the  initial  benchmark  location,  associated  ephemeris  data,  and 
output  the  position  of  each  spacecraft.  The  range  from  the  benchmark  to  the  spacecraft 
was  computed  and  the  residual  error  between  observation  and  measurement  time  delay 
produced.  The  tropospheric  propagation  delay  was  implemented  by  computing  residual 
errors  and  applying  a  geolocation  algorithm  to  generate  miss  distance  and  uncertainty. 

We  examined  the  relationship  among  benchmark  parameters  to  see  if  the  Choi 
model  had  a  predictably  favorable  impact  on  the  signal  geolocation.  Since  benchmark 
locations  are  known,  the  benchmark  miss  distance  is  a  powerful  indicator  of  the  system 
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model  accuracy.  We  would  expect  that  the  Choi  model  would  lead  to  a  reduced  miss 
distance  at  low  elevation  angle  signals.  The  next  parameters  of  interest  were  SMA  and 
SMI  uncertainty  which  reflect  the  statistical  uncertainty  in  the  position  solution.  Figures 
29  through  31  show  SMA,  SMI,  and  miss  distances  of  the  150  benchmark  signals  for  one 
channel  of  three  different  events.  These  figures  show  that  the  Choi  model  appears  to 
reduce  miss  distances  and  uncertainties  more  than  the  Hopfield  model.  In  Figures  29  and 
30,  one  can  find  that  some  SMA  parameters  (Choi  Model)  of  certain  benchmark  signals 
are  greater  than  those  of  the  Hopfield  SMA.  But  the  miss  distance  (Choi  model)  of  that 
benchmark  is  less  than  the  Hopfield  miss  distance.  These  cases  show  that  the  ellipse  error 
bound  using  the  Hopfield  model  never  covers  a  benchmark  location.  However,  the  ellipse 
error  boundary  using  the  Choi  model  is  greater  than  that  of  the  Hopfield,  always  covering 
the  benchmark  location.  Figure  32  shows  those  cases. 

Figure  33  shows  three  different  ellipse  errors  probability  (EEP  95%)  plots  for 
different  elevation  angle  ranges  in  one  benchmark  scan  data  set.  In  Figure  33.1,  the  ellipse 
indicates  the  case  of  the  elevation  angle  from  9.08  degrees  to  18.48  degrees.  The  miss 
distance  of  the  Choi  model  is  1,166.20  feet  and  the  Hopfield  model  3,570.85  feet.  In 
Figure  33.2,  the  ellipse  shows  the  case  of  the  elevation  angle  from  13.47  degrees  to  31.49 
degrees  and  then  decreased  to  15.66  degrees.  The  miss  distance  of  the  Choi  model  is 
355.17  feet  and  the  Hopfield  model  667.20  feet.  In  Figure  33.3,  the  ellipse  presents  the 
case  of  the  elevation  angle  from  17.77  degrees  to  negative  0.66  degrees.  The  miss 
distance  of  the  Choi  model  is  3,428.26  feet  and  theHopfield  model  4,406.05  feet.  Figure 
33.4  shows  six  ellipses  in  the  same  scale.  In  this  plot,  precise  geolocation  is  usually 
associated  with  the  scan  segment  received  at  the  high  elevation  angle  from  the  emitter  to 
the  spacecraft. 

The  geolocation  accuracy  of  the  Choi  tropospheric  model  is  improved 
approximately  50%  of  the  uncertainty  as  reflected  in  the  SMA  and  the  miss  distance.  We 
observed  that  the  Choi  tropospheric  time  delay  model  appears  to  improve  geolocation 
accuracy  of  the  benchmark  miss  distance  and  reduce  in  the  SMA  (uncertainty  in  the 
calculated  position).  We  can  infer  from  these  six  different  event  data  sets  that  the  Choi 
tropospheric  model  reflects  a  significant  improvement  in  the  benchmark  miss  distances  and 
uncertainties  (SMA  and  SMI). 


5.  Conclusion  and  Recommendation 

There  are  two  major  dominant  factors  in  tropospheric  RF  propagation  angular 
bending  and  time  delay.  These  measurement  errors  can  readily  be  corrected  by  means  of 
ray-bending  compensation  through  refractive  profiles  or  ray  tracing  algorithms  that  can  be 
computed  on  real  or  near  real  time  weather  data  such  as  surface  temperature,  pressure, 
and  relative  humidity  or  dew  point.  The  proposed  model  (Choi  tropospheric  model)  can 
generate  errors  of  both  range  and  time  delay  and  angle  of  arrival  in  near  real-time  and  real¬ 
time  for  calibration  of  tracking,  navigation,  and  satellite  location  in  the  world.  This 
enhancement  will  save  time  and  reduce  risk  in  locating  the  mission  target  and  search-and 
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Figure  29  SMA  for  one  channel  of  three  different  events 
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Figure  31  Miss  distance  for  one  channel  of  three  different  events 


Figure  32  Error  ellipse  plot  for  benchmark  CAIR 


rescue  (SAR)  mission  for  both  military  and  commercial  applications.  The  Choi  model  is 
currently  in  the  process  of  operational  test  and  evaluation  in  MRC  and  IAEP. 


The  Choi  tropospheric  propagation  model  clearly  affects  the  precision  geolocation 
process.  Geolocation  accuracy  of  benchmarks  is  improved  approximately  more  than  50% 
of  the  time  with  miss  distance,  SMA,  and  SMI.  Particularly,  we  would  expect  that  the 
Choi  tropospheric  model  would  lead  to  reduce  miss  distance  and  cover  the  benchmark 
location  in  the  lower  elevation  signals.  More  detailed  examination  of  the  scan  data  and 
emitter  geometry  may  help  improve  more  accuracy  of  geolocation  based  on  the  model. 

For  continuing  accuracy  improvement,  several  new  tasks  are  planned  to  implement 
a  total  atmospheric  model  (tropospheric  and  ionospheric).  The  following  tasks  will  be 
studied  in  FY2000. 


1)  Test  and  statistical  analysis  of  geolocation  accuracy  improvement  using  more  than 
1000  benchmarks  signal  data. 

2)  Examine  isochrone  correction  for  geolocation  using  TDOA  measurements. 

3)  Develop  a  general  atmospheric  model  by  combining  tropospheric  and  ionospheric 
propagation  characteristics. 

4)  Develop  a  frequency  dependency  characteristic  measurement  in  the  tropospheric 
region. 

5)  Develop  a  tropospheric  model  for  the  aircraft’s  angle  of  arrival  (AOA). 
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